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Abstract. The numerical hydrodynamic modelling of beat Cepheid behavior has been a long standing quest in 
which purely radiative models had failed consistently. We find that beat pulsations occur naturally when turbulent 
convection is accounted for in our hydrodynamics code. The development of a relaxation code and of a Floquet 
stability analysis greatly facilitates the search for and the analysis of beat Cepheid models. The conditions for 
the occurrence of beat behavior can be understood easily and at a fundamental level with the help of amplitude 
equations. 
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1. Introduction 

Double-mode or beat pulsations are not uncommon among 
RR Lyrae (the RRd stars) and classical Cepheids (c/. e.g. 
Kovacs 2001). Yet the numerical modelling of this phe- 
nomenon remained an unmet challenge for years. Partial 
success was achieved with purely radiative codes (Kovacs 
& Buchler 1993), but it took the incorporation of a model 
for turbulent convection into the hydrocodes to obtain 
results that are both robust and have astrophysical pa- 
rameters that are in agreement with the observations. 
Essentially simultaneously, but independently, KoUath et 
al. (1998) could model double-mode pulsations in Cepheid 
models and Feuchtinger (1998) in RR Lyrae models. In 
this paper we present an extension of this earlier work. 
We examine in some detail specific examples of both RR 
Lyrae and Cepheid models, as well as sequences of such 
models (in which L and M are held fixed and Toff is varied, 
thus mimicking approximately the horizontal branch and 
the Cepheid loops, respectively). 

2. Turbulent Convection Model Equations 

The turbulent convection model equations that we have 
incorporated in our hydrodynamics code are essentially 
those of Gehmeyr & Winkler (1992) and Kuhfufi (1986). 
They differ in the expression for the convective flux Fc 
from those of Stellingwerf (1982) and Bono & Stellingwerf 
(1994) and those of original paper (Yecko, KoUath & 
Buchler 1998), but are essentially the same as used by 
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Wuchterl & Feuchtinger (1998). To be specific we repro- 
duce here the equations that we use in our hydrodynamics 
code. 

The fluid dynamics part of the model calculations are 
given by the following equations: 
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The turbulent motion of the gas and the convection inter- 
acts with the hydrodynamics of the radial motion through 
the convective flux Fc, the viscous eddy pressure Pi, the 
turbulent pressure pt , and finally, through an energy cou- 
pling term C. The turbulent energy et is determined by a 
time dependent diffusion equation 
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The coupling term that connects the gas and the turbulent 
energy equations is given by: 



C ^ ad -J- [St - et) , 



(4) 



where A aA Hp, Hp = pr'^ /{pGM) is the pressure scale 
height, aA is the mixing length parameter. The a^'s are 
dimensionless parameters of order unity. Both the convec- 
tive flux and the source term of the turbulent energy (St) 
depend on the dimensionless entropy gradient 
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3.1. Time Dependence of Amplitudes and Phases 
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The remaining quantities are defined as 
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The turbulent convective equations can be derived on di- 
mensional and physical grounds with a number of dimen- 
sionless parameters (a's) that are of order unity but for 
which theory provides no numerical values. Our approach 
has been to attempt to calibrate them through a compar- 
ison of our results with the observational data. 

The calibration of the a's is a daunting task because of 
the large number of parameters on the one hand, and the 
large number of observational constraints, many of which 
require extensive full amplitude hydrodynamical model 
calculations. For example, important constraints are the 
locations and widths of the fundamental and first over- 
tone instability strips, the Fourier decomposition coeffi- 
cients of the light curve and of the radial velocity data, 
with their dependence on mctallicity. Such a global cali- 
bration has not yet been performed, and it is not even sure 
that it possible given the rather simple turbulent convec- 
tion recipe that we use. In that sense our results must be 
regarded as tentative. We stress however, that the types of 
behavior that we describe here are not sensitive to these 
a's, but that the masses, luminosities and Toff's at which 
they occur depend on them. In KoUath & Buchler (2001) 
we have presented typical modal selection diagrams for 
Cepheid and RR Lyrae models (figures 14 and 15) that 
show the calculated pulsational behavior of the models in 
a Hertzsprung-Russell diagram. 

3. Transient Behavior of Models 

The first issue that we address concerns the confirmation 
that a given model is indeed undergoing stable double- 
mode pulsations rather than being involved in a switch 
from one mode to another. Deciding on the basis of a sin- 
gle hydrodynamical calculation whether a model is under- 
going steady double-mode pulsations can be fraught with 
peril. The model may give the false impression of having 
achieved steady behavior when it is actually still in a tran- 
sient state even after thousands of periods of integration, 
and it will end up pulsating with a single frequency. The 
reason for this will become clear shortly. 



One of the several tools that are needed is the calculation 
of the time dependence of the amplitudes and phases of 
the modes that are excited in the models. Instead of us- 
ing a time dependent Fourier analysis as advocated by 
Kovacs, Buchler & Davis (1987, KBD87) we now use 
the much more powerful and accurate analytical signal 
method of Gabor (1946) to reconstruct the time depen- 
dent amplitudes and frequencies of the pulsation modes, 
(c/. e.g., Cohen 1994). The method can be extended to 
multi-component signals with the help of filters that re- 
strict the power to the desired frequency components. It 
is convenient to make this filtering in Fourier space where 
it involves just a product, and to combine it with the def- 
inition of the analytical signal (KoUath & Buchler 2001) 



Zk{t) = afc(t)e^'^'=(*) 
1 



dt's{t')e~ 
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where s{t) is the real valued signal that is to be analyzed, 
H{u} — ujk) is the filtering window that is centered on the 
desired Uk- A Gaussian window with a half width of Wfc/lO 
generally gives satisfactory results for the stellar pulsa- 
tions. The resulting amplitudes and phases describe then 
the transient behavior of the given mode. The analytical 
signal approach yields a much better resolution than the 
KBD87 approach because it is not necessary to average 
the amplitudes in time to get smooth results. 

With the analytical signal method it is therefore very 
easy and fast to define unambiguously the instantaneous 
phase ip{t) and amplitude A{t) of a signal. 

3.2. Phase Plots 

The computed time dependent amplitudes, in our case 
A(){t) and Ai{t), can then be used to visualize the tran- 
sient behavior of the pulsating model in an {Aq,Ai) phase 
space. Our amplitudes are defined as those of the nor- 
malized radial displacement of the photosphere, 6Rph/R*, 
where i?, is the photospheric radius of the equilibrium 
model. 

In Fig. on the left we display the results of a num- 
ber of hydrodynamical calculations made for the same 
Galactic Cepheid model (M=4.7 Mq, L=1335 Lq, T^s 
= 5580 K, Z=0.02), but started (kicked from the unstable 
equilibrium) with different initial conditions. The tracks 
all start with small amplitudes (lower left side). 

Several notable features stand out in the left of Fig. |[ 
Since the thick dots denote points that are equally sep- 
arated in time one sees that the models all evolve very 

^ Physicists are generally familiar with this analytic sig- 
nal concept through the Kramers-Kronig dispersion relations 
(Jackson 1975). If s{t) represents the real part of an assumed 
complex analytical function Z{t), then the imaginary part s{t) 
of Z{t) can be obtained via a Cauchy integral, which through 
contour deformation becomes a Hilbert transform, which in 
turn can be converted into a one-sided Fourier transform. 
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Fig. 1. Cepheid model with T^s^ddSOK. 

Left: The flow in the (Aq, Ai) phase space. The large dots represent the hydrodynamical results at equal time intervals, 
so that the spacing of the dots gives an indication of the speed. The short lines represents the flow field (dots denote 
the base of the normalized vectors). 

Right: The stable/unstable fixed points are denoted by flUed/open circles. The integral lines (heteroclinic connections) 
that connect the fixed points are shown as solid lines with arrows indicating the evolution direction. The thick line 
is also the separatrix, i.e. the boundary between the basins of attraction of the two stable fixed points, viz.the stable 
DM and the fundamental. 



quickly toward an arc along which the transient evolution 
is extremely slow. If only one model were computed this 
'convergence' could easily be misinterpreted as a steady 
double-mode pulsation. The figure indicates that the point 
on the ^0 axis with {Aq = 0.026, Ai = 0) is an at- 
tractor. The three rightmost tracks transit to this funda- 
mental attractor after a sufficiently long integration time. 
This corresponds of course to the fundamental limit cy- 
cle. Similarly the point on the arc with {Aq = 0.023, Ai = 
0.007) repels the tracks, and it corresponds to an unstable 
double-mode pulsation. Further up along the arc a second 
attractor with (^o = 0.015, Ai = 0.014) is clearly visible. 
Six of the tracks converge toward it from the left and three 
of them from the right. This attractor corresponds to sta- 
ble double-mode pulsations. Finally, one notes that the 
single-mode overtone limit cycle at (^o — 0,Ai = 0.017) 
repels because it is unstable. All tracks run away from the 
origin because the equilibrium model is linearly unstable 
in both modes. 

We shall see that this whole picture can be clearly un- 
derstood after we have introduced amplitude equations. 
They allow us to draw the flowfield that has been super- 
posed on the left-hand figure. 

3.3. Amplitude Equations 

It has been known for a while now that one can de- 
velop a global understanding of the 'modal selection' prob- 



lem with the help of the amplitude equation formalism 
(Buchler & Goupil 1984). The beauty and power of the 
amplitude equations are their simplicity and their very 
basic and generic nature. They apply to any system, be it 
a pulsating star, a biological system, or any other dynam- 
ical system, in which (a) the relative growth-rates of the 
excited modes are small; (b) the pulsations are sufficiently 
weakly nonlinear so that an expansion in the amplitudes 
is justified. Only selected powers of the amplitudes appear 
(Buchler & Goupil 1984, Buchler 1993), and the omitted 
terms are unessential for the modal selection (bifurcation 
diagram), i.e., the nature of the dynamical behavior does 
not depend on them. 

Condition (a) is clearly satisfied for RR Lyrae and clas- 
sical Cepheids, since typically the growth-rates are of the 
order of a percent of the frequencies. That condition (b) is 
also satisfied manifests itself, for example, by the fact that 
the nonlinear periods are very close to the linear ones. 

The form that the amplitude equations take on de- 
pends on whether or not there are nearby resonances. ^ 
For the beat RR Lyrae there is no low-order resonance in 
the period range of interest, and for the beat Cepheids the 
nearest resonance is too far to be influential. The appro- 
priate amplitude equations for the nonresonant case are 

do = ao (do + Qoo|aoP + QoilaiP + •^'olaoPlaiP 

^ By resonance we mean a relation of the form niujo + n2UJi ~ 
or niijJo+n2LOi+n2ijJ3 ~ 0, where ni,n2, na are small positive 
or negative integers. 
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Fig. 2. Left: Behavior of the amplitudes of single-mode (SM) and double-mode (DM) pulsations along a constant 
L convective Cepheid model sequence; Solid lines: stable pulsations, dashed lines: unstable pulsations; c/. text. The 
filled and open circles locate calculated hydro-models. 
Right: RR Lyrae sequence. (Amplitude in units of SRph/R*-) 




+ i?oo|ao|^ + ^01 (14) where the dot represents the time-derivative and aj are the 

linear eigenvalues for an assumed exp(crt) dependence, 

aj = Kj + iuj (16) 

^ ^ 2 2 Ti^G flj 's are the complex amplitudes of the two excited 

oi = ai (cti +Qio|ao| +Qii|ai| +S'i|ao| |ai| modes, and Qjk, S-j and Rju are the complex nonlinear 

+ i?io|ao|^ + i?ii|ai|'') (15) couphng constants. 
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It is found to be more convenient to use real am- 
plitudes Aj{t) and phases ^j(t), defined by aj{t) = 
Aj{t) e'}q)iipj{t), because this produces a complete decou- 
pling of the amplitudes from the phases. One finds 



+s,AlAl + r,oAt + r,,At) (17) 

and 

ifiO ^ UJO + qooAl + qoiAf + SqAIAI (18) 

ifii = uji + gioA^ + qiiAj + siAlAl (19) 
where 

Qjk = qjk +iqjk (20) 

Sj = Sj +isj (21) 

R]k = Tjk +irjk- (22) 



We have disregarded the imaginary parts of the higher 
order terms in the phases because they are unimportant. 
One also notes that all equations can be expressed in terms 
of the squares Aq and Af . 

The constant amplitude solutions of Eqs. (|lj, [l^) the 
so-called fixed points are therefore obtained from solving 
the equations 



Ko + qooAl + qoiAj + soAlAj 

+rooA^o + roiAt ^ (23) 

^1=0 (24) 

Ki + qwAl + qiiAl + siAlAl 

+rio4 + siiAf = (25) 

^0 = (26) 



pairwise between the two sets. They correspond to steady 
nonlinear pulsations. In the earliest use of amplitude equa- 
tions for describing DM pulsations (Buchler & Kovacs 
1986; KoUath et al. 1998) the amplitude equations were 
truncated at the lowest nontrivial, i.e. the cubic terms. It 
was shown that in this order the single-mode fixed points 
and steady DM pulsations cannot simultaneously be sta- 
ble for the same stellar model. Since then it has been re- 
alized that the behavior of the hydrodynamical models 
and the observational constraints imposed by the Beat 
Cepheids and RR Lyrae stars can be a little more com- 
plicated than the cubic nonresonant amplitude equations 
allow. However, the more complicated behavior can be 
captured if we keep the next order nontrivial terms, viz. 
the quintic ones in the amplitudes (Buchler et al. 1999). 
It is not necessary though to keep all quintic terms, and 
for simplicity, these authors retained only the roo and rn 
cubic terms and discussed the fixed points of the ampli- 
tude equations. Here we keep instead the sq and Si cubic 
terms because they give a slightly better description of the 
results of the hydrodynamical calculations. 



3.4. Determination of the Nonlinear Coefficients 

We can use the hydrodynamical calculations that are re- 
ported in Fig. |l| to determine the values of the unknown 
nonlinear coupling coefficients that appear in the ampli- 
tude equations. This is done as follows: The analytical 
signal method gives the time dependent amplitudes suffi- 
ciently accurately so that the derivatives A can be used di- 
rectly in Eqs. |l^ to obtain the linear growth rates and the 
nonlinear coupling coefficients in these equations through 
a linear least squares fit. (With the method described in 
Buchler & Kovacs (1987) it was necessary to make a te- 
dious nonlinear fit.) 

Once the coefiicients are known, we can then produce 
a vector plot u = {Aq,Ai). The computed vector plot 
has been superposed on the transient evolutionary tracks 
in Fig. |l|, and gives a global view of the flow field. (The 
dots represent the bases of the vectors) . The speeds have 
been normalized because they span a very broad range of 
values. 

The fixed point scenario is summarized in the right side 
of Fig. 1^. The solid circles denote the stable fixed points 
and the open circles the unstable ones. We recall that the 
five fixed points each correspond to a steady pulsational 
behavior. The origin represents the unstable equilibrium 
model (an unstable node point; for a mathematical defini- 
tion of these terms c/. Lefshetz 1977), the two fixed points 
on the axes are the stable fundamental and unstable first 
overtone limit cycles (a stable node point and a saddle 
node point, respectively). The remaining two fixed points 
are the stable and unstable double-mode pulsations (again 
a stable node point and a saddle node point, respectively). 
Also shown in the figure are the integral lines that link the 
fixed points (the heteroclinic connections), and the special 
one, the separatrix (thick line), that separates the two at- 
tractor, i.e. the regions of initial conditions that lead to 
the stable DM pulsation and to the fundamental limit cy- 
cle, respectively. 

The amplitude equations give an accurate and global 
picture of the pulsational behavior of the models. 

4. Sequences of Models 

4.1. Turbulent convective models 

In this section we examine constant L and M sequences 
of turbulent convective models, one for Cepheids and one 
for RR Lyrae stars. These sequences display a somewhat 
different behavior, but that has more to do with the chosen 
stellar parameters L and M than with the nature of the 
stars. In fact both Cepheids and RR Lyrae exhibit both 
types of behavior in some range of L and M. 

4.1.1. Cepheid Sequence: 

In Fig. H we report the results of the convective Cepheid 
model sequence (M=4.7 Mq, L=1335 Lq) with Toff vary- 
ing from 5690 to 5530 K. The solid circles correspond to 
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TABLE 1.- Nonlinear Coupling CoefRcients and Growth Rates. 



Convective Cepheid Model sequence. 



5560 


-0.92 


-1.49 


-5.77 


-11.99 


-3436.0 


2449.0 


0.00066 


0.00329 


5580 


-0.93 


-1.50 


-5.62 


-12.17 


-2816.0 


4994.0 


0.00065 


0.00351 


5600 


-0.94 


-1.45 


-5.48 


-12.03 


-2653.0 


6594.0 


0.00063 


0.00370 


5640 


-0.98 


-1.63 


-5.14 


-11.53 


-2054.0 


5567.0 


0.00060 


0.00408 


5680 


-1.03 


-1.87 


-4.66 


-10.91 


-1545.0 


830.7 


0.00057 


0.00441 



Convective RR Lyrae Model sequence. 



6480 


-3.08 


-10.28 


-19.81 


-44.26 


-2701.0 


-1000.0 


0.00964 


0.03953 


6500 


-3.09 


-10.18 


-20.24 


-43.98 


-2826.0 


-545.7 


0.00936 


0.03961 


6505 


-3.07 


-10.13 


-20.35 


-43.94 


-2750.0 


-500.4 


0.00924 


0.03967 


6515 


-3.08 


-10.08 


-20.62 


-43.90 


-2833.0 


-383.7 


0.00910 


0.03976 


6525 


-3.07 


-9.99 


-20.75 


-43.72 


-2792.0 


-359.9 


0.00890 


0.03971 


6600 


-3.22 


-9.92 


-22.26 


-42.38 


-2781.0 


36.2 


0.00774 


0.03866 


6700 


-3.85 


-10.08 


-22.29 


-40.16 


-2012.0 


-3303.0 


0.00605 


0.03519 



the amplitudes of the computed stable DM pulsators, and 
the open circles to single-mode first overtone or funda- 
mental pulsators. 

In Table la we show the nonlinear coupling coefficients 
for the Cepheid sequence. The calculations are laborious 
albeit straightforward, because for each model in the se- 
quence to the numerical hydrodynamics calculations have 
to be made for a set of different initial conditions, and then 
the coefficients have to be computed through a nonlinear 
fit. One notes that the coefficients do not vary much from 
one model to another. This, by the way, serves also as a 
confirmation of the quality of our fits. 

The lines in Fig. || show the amplitudes of the possible 
pulsational states of the models. The lines were obtained 
by computing the fixed points of the amplitude equations 
with average values of the nonlinear coupling constants, 
but using the TeS dependent growth rates kq and ki. 

The solid lines denote stable and the dashed lines un- 
stable behavior. The first overtone limit cycle is seen to 
be stable down to T^s = 5647 K. The fundamental limit 
cycle is stable up to Tcs— 5598 K, and a stable DM pul- 
sation exists between T^s = 5647 and 5568 K. (A second 
DM pulsation is also possible between Tes = 5598 and 
5568 K, but is unstable and therefore astronomically not 
achievable.) The sequence thus shows hysteresis, i.e. in the 
range Toff = 5598 to 5568 K the pulsational state depends 
on whether the stellar evolution is blueward or redward. 

Because Eqs. (|2^ ^^6|) are a function of the squared 
amplitudes, another visualization of the bifurcation sce- 
nario is afforded by Fig. || which shows the fixed points 
in an {Aq,A1) plot for the same Cepheid model sequence. 
The loci defined by Eqs. ( p3| , [25| ) are shown as dashed and 
solid lines, respectively. The remaining two loci are the 
positive X and y axes. The solid and open circles represent 



stable and unstable fixed points which are the intersec- 
tions of these loci. 

4.1.2. RR Lyrae Sequence: 

The results of this sequence are reported in Fig. ^. Here we 
have the possibility of a fundamental limit-cycle with Tcs 
< 6700K, of an overtone limit-cycle with Teff > 6500 K. 
The Teff range of DM behavior is extremely narrow. Again 
we have hysteresis. Along a redward track the model starts 
off as an RRc star, then becomes a beat RR Lyrae (RRd) 
and finally changes into an RRab. Along a blueward ex- 
cursion the model starts of a RRab and then turns into 
an RRc, without ever being RRd. 

Compared to Fig. ^ the situation is now reversed in 
terms of the relative locations of points a and 7. Again, 
if we could disregard the quintic terms for the sake of 
argument, points a and /3 would merge, and in the sce- 
nario described by Buchler & Kovacs (1986), we are now 
in an 'either fundamental or overtone' situation (X* < 0), 
as confirmed by the values in Table lb for a sequence of 
RR Lyrae (M=0.77 Mq, L=50.0 Lq, Z=0.0001) models. 
The DM behavior occurs here only because of the looping 
back due to the quintic terms, and as a result is much 
narrower than in the Cepheid case. 

Fig. ^ shows the fixed points in an {Al,Al) plot for 
the same RR Lyrae model sequence. The loci defined by 
Eqs. (^, ^5|) are shown as dashed and solid lines, respec- 
tively. The remaining two loci are the positive x and y 
axes. The solid and open circles represent stable and un- 
stable fixed points which are the intersections of these loci. 

We summarize the modal selection results obtained for 
a number of RR Lyrae sequences in Figure ^ which shows 
the regions in which the various pulsational behavior oc- 
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Fig. 4. Modal selection of RR Lyrae model sequences The gray shading of the modal states are described at the top of the left 
figure: F: fundamental mode, F/DM either fundamental or double mode, DM double mode, F/01 either fundamental or first 
overtone, 01 first overtone 



curs. In these calculations the metallicity was at Z=0.0001, 
since according to our calculations the metallicity has only 
a minor effect of the RR Lyrae edges. All blue and red 
edges are nonlinear ones except for the F red edge, which 
is from linear stability analysis. Note that at low lumi- 
nosity and high mass an additional F region appears. We 
caution again that these results depend on our choice of 
turbulent convective parameters and that they may have 
to be revised, at least quantitatively, when a more defini- 
tive calibration is performed. 

5. Discussion 

We have already pointed out that radiative models have 
failed to yield steady DM pulsations, at least with a behav- 
ior that satisfies the observations. It is therefore of interest 
to dwell a little on why convective models do. 

First we note that the looping back between 5600 
and 5570K in the Cepheid models and between 6510 and 
6500K for RR Lyrae models in Fig. || occurs because of 
the presence of quintic terms in the amplitude equations. 
If for the sake of argument we were to ignore the quintic 
terms then the a and 7 points in Fig. ^ would merge, and 
would be to the right of the point a, where the fundamen- 
tal amplitude vanishes. 

In Figs. I the curves would turn into straight lines, as 
shown schematically in Fig. |[ In the situation depicted on 
the left side the model pulsates in either the F or the 01 
limit cycle, depending on the past evolutionary history. 
The right side corresponds to the DM case. It is easily 
shown (Buchler & Kovacs 1986) that a necessary condition 
for double-mode behavior is that V = qooqu — goi^io > 
0. This condition is never found to be satisfied in purely 
radiative models in which the cross-coupling terms qoi and 
qio always dominate over the self saturation terms qoo and 

One notes in Fig. |that the rcff=5610K Cepheid model 
has the same topology as the right figure with the DM 
behavior, and that the Teff=6530K RR Lyrae model has 
that of the left side. Table 1 shows that indeed V > Q for 




Fig. 5. Schematic plot. Solid circles: stable fixed points, 
open circles: unstable fixed points. Left: F or 01 regime, 
2? < and Right: DM regime, V > 0. 

the Cepheid model sequence and V < for the RR Lyrae 
one, even though this criterion applies only approximately 
when the quintic terms arc considered. Nevertheless it ex- 
plains why the DM region in the convective Cepheid model 
sequence is much broader than that of the RR Lyrae (it 
would exist even in the absence of the quintic terms) . 

In KoUath et al. (1998) we had somewhat hastily sug- 
gested that the abrupt disappearance of the DM solution 
was probably due to a pole caused by the vanishing of "D. 
There is indeed a pole in the broad vicinity, but it is too 
far to play a role. Instead, as we have shown here, the be- 
havior is due to the appearance and collision of two DM 
fixed points, a situation that is made possible because of 
the importance of the quintic terms which generate the 
two DM fixed points. 

6. Conclusions 

We have shown that there are two reasons for the oc- 
currence of DM behavior with the turbulent convective 
models. One is the increase of the cubic saturation terms, 
goo and qu, relative to the cross-coupling terms, goi and 
(jiQ. This is a rather subtle and hidden effect of turbulent 
dissipation because these coefficients are given as compli- 
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cated integrals over the structure of the star (Buchler & 
Goupil 1984). 

The other reason is the necessity to include quintic 
terms in the description which, in turn, enrich the com- 
plexity of the bifurcation diagram, in particular by caus- 
ing the simultaneous occurrence of two DM fixed points - 
compare Figs. |^ and ^. 

One of the interesting questions that we have not ad- 
dressed concerns the time it takes to switch from one pul- 
sational mode (fixed point) to another one. This problem, 
in our opinion, had never been correctly addressed in the 
pulsation literature, and we discuss it in a companion pa- 
per (Buchler & KoUath 2001) where we show that this 
transition does generally not occur on the thermal time 
scale tth — K~^, but on a much longer time scale, which, 
depending on the type of bifurcation is either set by the 
stellar evolution time scale, more precisely {dn/ dt)~'^'^ ^ or 
by the thermal time tth, however generally multiplied by 
several orders of magnitude. 
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